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Abstract. We present morphological and spectroscopic studies of the northwest rim of the supernova remnant 
RX J1713. 7-3946 based on observations by the Chandra X-ray observatory. We found a complex network of nonthermal 
(synchrotron) X-ray filaments, as well as a 'void' type structure - a dim region of a circular shape - in the northwest rim. It 
is remarkable that despite distinct brightness variations, the X-ray spectra everywhere in this region can be well fitted with a 
power-law model with photon index ranging T = 2.1-2.5. We briefly discuss some implications of these results and argue that 
the resolved X-ray features in the northwest rim may challenge the perceptions of the standard diffusive shock-acceleration 
models concerning the production, propagation and radiation of relativistic particles in supernova remnants. 
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1. Introduction 

The supernova remnant (SNR) RX J 17 13.7-3946 has proven 
to be a prominent source of nonthermal X-rays and presum- 
ably also y-rays, thus providing strong evidence that shell- 
type SNRs are sites of acceleration of galactic cosmic rays. 
This source was discovered during the ROSAT All-Sky Survey 
(Pfeffermann & Aschenbach ll996> . Observations with ASCA 
have revealed intense synchrotron X-ray emission from the en- 
tire remnant (Koyama et al. 119971 Slane et al. I1999> . Unlike 
SN 1006, no evidence for thermal X-ray components has yet 
been found. 

At the north perimeter of RX J1713. 7-3946, a molecu- 
lar cloud (cloud A) has been found with a high CO (7=2- 
l)/CO (7=1-0) ratio that suggests possible interaction between 
the cloud and the SNR shell (Butt et al. l2"0TTTY The cloud has 
positional association with the unidentified y-ray source 3EG 
J1714-3857 (Butt et al. l2001i . Most recently, Uchiyama et al. 



( 2002 1 reported an unusually flat-spectrum X-ray source (AX 
J1714. 1-3912) coincident with this cloud, and argued that the 
flat spectrum can be best interpreted by bremsstrahlung from 
either sub-relativistic protons or mildly relativistic electrons. 
The high-energy particles responsible for the X-ray and y-ray 
emission from this cloud are likely associated (in one way or 



another) with the SNR-cloud interaction. For the preferred dis- 
tance to SNR RX J1713.7-3946 of d 6 kpc (Slane et al. 
1999 1, based on the kinematic distance to cloud A, the age of 
the remnant is estimated to be > 10000 yr. A younger age of 
~ 2000 yr cannot be, however, excluded if d = 1-2 kpc. 

The CANGAROO collaboration (Muraishi et al. 2000) re- 
ported the detection of TeV y-ray emission from the direction 
of the northwest (NW) rim, the brightest region of synchrotron 
X-rays in RX J1713. 7-3946. If confirmed, the TeV y-radiation 
would provide direct and unambiguous evidence for the pres- 
ence of particles (electrons and/or protons) accelerated to very 
high energies. Recently Enomoto et al. (2002i published the 
spectrum of TeV emission based on the new CANGAROO ob- 
servations. The spectrum is claimed to be quite steep with a 
power-law photon index F = 2.8 + 0.2 between 400 GeV and 
8 TeV. They argued that the steep spectrum is inconsistent with 
the inverse Compton (IC) model, but could be explained by n°- 
decay gamma-rays. If true, for the canonical shock-acceleration 
spectrum of protons with power-law index s ~ 2, this would 
imply a cutoff energy in the proton spectrum around 100 TeV. 
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Subsequently, Reimer & Pohl (2002i and Butt et al. (2002i ar- 
gued that this interpretation would violate the y-ray flux up- 
per limits set by the EGRET instrument. However, this is not a 
sufficiently robust argument to be used to dismiss the hadronic 
origin of the reported TeV emission. Adopting a slightly harder 
proton spectrum, e.g. with spectral index s < 1.9, it is possible 
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to avoid the conflict with the EGRET data. Even for proton 
spectra steeper than s = 2, it is still possible to suppress the 
GeV y-ray flux, if one invokes the effects of energy-dependent 
propagation of protons while traveling from the accelerator 
(SNR shocks) to the nearby clouds (see e.g. Aharonian l2001> . 
Moreover, the lack of GeV y-rays can be naturally explained 
by confinement of low-energy (GeV) protons in the supernova 
shell, in contrast to the effective escape of high-energy (TeV) 
protons. 

On the other hand, the arguments against the IC model of 
TeV emission should be backed by thorough theoretical stud- 
ies based on higher quality data from the radio, X-ray and y-ray 
domains. X-ray observations are of particular interest because 
the synchrotron X-ray spectra reliably "control" the predictions 
of IC emission at TeV energies. In this paper, we report on 
the X-ray study of the NW rim of SNR RX J1713.7-3946 us- 
ing archival data obtained with the Chandra X-ray Observatory 
(Weisskopf et al. I1996i . We demonstrate that the X-ray emis- 
sion has remarkable substructure with bright filaments accom- 
panied by curious dark voids. The observed features set new 
standards in X-ray studies of SNRs which should help to under- 
stand deeper the nonthermal processes of particle acceleration, 
propagation and radiation in supernova shocks. 

2. X-ray data 

The observation of the northwest (NW) rim of the SNR 
RX J1713.7-3946 was carried out on 2000 July 25 (Obsld 
736) with the Advanced CCD Imaging Spectrometer (ACIS) on 
board Chandra. The data were obtained with the ACIS-I array 
consisting of four front-side-illuminated CCD chips. As shown 
in Fig.Q] the array was placed to cover the brightest portion of 
the remnant. We analyzed the processed level 2 event data (pro- 
cessing version R4CU5UPD14.1), retrieved from the Chandra 
public archive, with a new gain map applied. The overall light- 
curve was examined for possible background flares or fluctu- 
ations; no significant time-variation was found. During such 
a quiescent period, the instrumental background is relatively 
constant (+10 %) with time 1 . Indeed, the count rate in the 8.5- 
10 keV energy interval is typically only several percent lower 
than the (standard) blank-sky data sets. Since this difference 
is within the statistical uncertainties of our results presented 
in this paper, the instrumental background component is taken 
from the blank sky data. After data processing, we obtained an 
effective exposure time of 29.6 ks. 

2.1. Image analysis 

Figure |2] shows the Chandra image of the NW rim of SNR 
RX J1713. 7-3946 in the energy range of 1-5 keV with the 
photon counts summed in pixels of 4" x 4". Point sources are 
removed, based on a wavelet decomposition using the CIAO 
software wavdetect (Freeman et al. 2002 1. We extracted the 
profile of photon counts from the rectangular region depicted 
in Fig.|2l summed over the vertical narrow dimension of the 
rectangle (Fig.|5|l. Figure |4] shows the smoothed X-ray images 

1 For details, see http://cxc.harvard.edu/contrib/maxim/bg/ 
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Fig. 1. Overall ASCA view of the SNR RX J17 1 3.7-3 946 re - 
gion in the hard 5-10 keV energy band (Uchiyama et al. 2002 1, 
and the Chandra's 17' x 17' field-of-view. The SNR has an an- 
gular radius of about 40'. A transient X-ray source (Koyama et 
al. 1997 1 is blanked out. Depicted are the 68% confidence error 
contour of the EGRET unidentified source 3EG J17 14-3857 
(Butt et al. 2001 ), the approximate outline of the cloud A, and 
the center of the cloud C (a triangle). 
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Fig. 2. Chandra ACIS-I image of the NW rim of 
RX J1713.7-3946 in the broad energy band (1-5 keV) 
with the photon counts (linear scale from to 15 counts) 
accumulated in 8 x 8 pixels. The coordinates are R.A. and decl. 
(J2000). 



in the soft 1-3 keV and hard 3-5 keV bands. The images have 
been adaptively smoothed with the CIAO software csmooth, 
where the local photon counts are utilized to determine the 
width of the Gaussian smoothing kernel at each position. The 
map of kernel widths was calculated based on the broad 1-5 
keV image. In Fig. [5] we show the spatial distribution of the 
hardness ratio, (3-5 keV)/(l-3 keV), of the smoothed images, 
after corrections for the mirror vignetting and exposure time. 

Dozens of compact and extended emission features have 
been detected throughout this region. The ASCA image of the 
NW rim now is resolved to distinctive components. Chains of 
bright filaments are evident within a more diffuse and fainter 
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Fig. 3. Horizontal profile of the photon counts (1-5 keV) for 
the rectangle region shown in Fig.|2] integrated over the vertical 
dimension. Distance is in units of 4", from east to west. 
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Fig. 4. Smoothed image of the NW rim of RX J1713.7-3946 
in a) soft 1-3 keV band and b) hard 3-5 keV band, obtained 
with four (each with 8.'4 x 8'4) CCD chips. The scales are in 
units of counts s _1 arcmkr 2 . The corrections for the inter-chip 
gaps and vignetting have not been applied. 




Fig. 5. The spatial distribution of the hardness ratio, (3-5 
keV)/(l-3 keV), in the NW rim, obtained by dividing the 
smoothed brightness map (in units of photons s _1 arcmin" 2 ) in 
the 3-5 keV energy band to that in the 1-3 keV band. 




Fig. 6. The regions selected for the spectral analysis, overlaying 
the same image as Fig.|2 The rectangular regions were selected 
for the bright filaments (1^4), plateau (5-7), and dark voids (8- 
10). 



plateau. The surface brightness of the filaments is on average 
about two times higher than that of the surrounding plateau, 
although the latter constitutes the bulk (~ 80 %) of the X-ray 
flux detected from the rim region comprised of the filaments 
and plateau. The most prominent and largest filament located 
near the center of the image and running from north to south, 
has an apparent width of 20", which corresponds to linear 
scale of AR ^ 0.58 d& pc, where d(, is the distance in units of 6 
kpc. Adjacent to the west side of the prominent filament, there 
can be seen an interesting morphology - a "void", i.e. a dim 
circular region with a radius of 3'. South of this circular void, 
we can see another void structure which seems to be extending 
beyond the observed field. 

The features mentioned above are observed both in the soft 
and hard images. The two images are strikingly similar, thus 
implying small spectral variations across the NW rim including 
the filament, plateau, and dark void regions. This is quantita- 
tively shown in the hardness-ratio map (Fig.|5}- The northwest 
corner of the observed field, outside the boundary of this SNR, 
shows relatively hard emission, which would originate in hot 
optically-thin thermal plasmas pervading the Galactic plane. 

2.2. Spectral analysis 

We have examined the X-ray spectra extracted from various re- 
gions in the NW rim, including the filaments, plateau, and void 
regions (Fig.|6jl. The instrumental background components are 
subtracted by using blank-sky observations. Since the contri- 
bution of diffuse Galactic X-ray emission to the spectra is es- 
timated to be about 5 % and 10 % for the filament and plateau 
regions, respectively, it has negligible effects on the spectral 
properties derived from these regions. Note, however, that for 
the void regions, the contamination could go up to 20 %. 

Like the spatially integrated spectrum that has been discov- 
ered by ASCA, the spectra from individual regions can be rep- 
resented as featureless continuum. This can be best interpreted 
by synchrotron emission produced by relativistic electrons with 
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Fig. 7. ACIS-I spectrum of the filament region 2. Superposed 
on the data points is histogram of the best-fit power-law model 
with interstellar absorption. 
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Fig. 8. Results of spectral fits with an absorbed power-law 
model for several selected regions: the bright filaments (region 
1-4-), plateau (region 5-7), and dark voids (region 8-10). From 
top to bottom panels, plotted are the best-fit values (with their 
90% errors) of photon index F and absorbing column den- 
sity A^h in units of 10 22 crrT 2 , surface brightness in units of 
10" 12 erg cirT 2 s" 1 arcmin" 2 , and reduced chi-square^- 2 . 



multi-TeV energies. We fit the spectrum in the 0.8-7 keV in- 
terval individually with a power-law model taking account of 
photoelectric absorption along the line-of-sight, and obtain sta- 
tistically acceptable fits in all cases. In Fig.0 we present the en- 
ergy spectrum extracted from region 2 as a characteristic exam- 
ple, together with the best-fit model. We show the spectral pa- 
rameters for these regions in Fig. [8] It is demonstrated that the 
best-fit values of photon index T [F(s) = Ke r (photons cm" 2 
s _1 keV" 1 )] and absorbing column density A^h are remarkably 
similar to each other, despite substantial brightness variations, 
with average values of F ^ 2.3 and A^h - 0.9 x 10 22 cm" 2 . 
These parameters agree fairly well with the reported ASCA 
spectrum integrated over the NW rim: F = 2.4 ±0.1 and 
N H (0.81 + 0.06) x 10 22 cm" 2 (Koyama et al. fT997l . 



Generally, a simple power-law function provides very con- 
venient and meaningful parametrizations, especially for non- 
thermal X-ray in limited energy intervals. However, the syn- 
chrotron origin of X-ray emission a priori implies the presence 
of spectral cutoffs associated with the unavoidable (see below) 
cutoff in the acceleration spectrum of parent electrons. More 
specifically, in the case of SNRs the spectrum of synchrotron 
X-rays is formed by the highest-energy electrons distributed 
in the cutoff region. Although, the current models of diffusive 
shock acceleration do not provide definite predictions about the 
spectral form of electron distribution in the cutoff region, an ex- 
ponential cutoff is commonly believed to be a natural assump- 
tion. Below we will adopt a more relaxed (generalized) form 
for the electron distribution, namely 



N(E) oc E~ s exp 



(1) 



with /3 =1/2, 1, and 2, where E m is the maximum energy of 
accelerated electrons. 

For this energy distribution of electrons, the ^-functional 
approximation gives a simple analytical presentation for the 
differential spectrum of synchrotron radiation: 



F(s) oc s r exp 



(2) 



where a = /3/2, r = (s + l)/2 is the photon index and bq is the 
cutoff energy. Generally, so = fe c , where 



5.3 M- H-^M keV 
UOjuG/llOOTeV/ 



(3) 



is the characteristic energy of synchrotron photons emitted by 
an electron of energy E = E m , and £ is a parameter introduced 
to adapt this presentation to the results of the accurate numeri- 
cal calculations. These calculations show that for ( = 1 (but not 
for £ a 0.3 as often used in literature), Eq. @ gives an adequate 
description of the synchrotron spectrum with accuracy < 25 % 
in a rather broad energy band, including the cutoff region (e.g. 
Aharonian 2000). 

Therefore below we invoke the spectral presentation given 
by Eq. @ in order to estimate the allowed range for the cutoff- 
energy position in the Chandra data. We treated the cutoff en- 
ergy eo as a free parameter, and derived the best-fit values for 
different combinations of fixed parameters F and a, namely 
r = 1.5, 1.75,2.0 and a = 1, 1/2, 1/4. We note that for the 
standard shock-acceleration spectrum with s — 2, one should 
expect F = 1.5, provided that the spectrum of electrons is not 
suffered deformation due to radiative and non-radiative (e.g. 
escape) losses. However, typically the cooling time of highest 
energy electrons responsible for synchrotron X-ray emission is 
shorter than the source age, even for young, 1000 yr old SNRs. 
This leads to the spectral steepening, thus, F ~ 2.0 (correspond- 
ing to s — 3) would be a more realistic choice. 

We have performed spectral fits for one of the filamentary 
regions, region 2, by the model defined by Eq. (0 together with 
interstellar absorption. We have obtained good fits (xl 1 .0, 
y = 127) for each set of parameters, the absorbing column den- 
sities being found in the range A^h = (0.66-0.79) x 10 22 cm" 2 . 
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Table 1. Cutoff photon energy obtained for the synchrotron fil- 
ament (region 2). 
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5 2 +L1 
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1/2 








a = 


1/4 


0.08^ 


r = 


1.75 


a = 


1 


8 4 +34 






a = 


1/2 


5 +56 






a = 


1/4 


44+ 18 

"•^-0.26 


r = 


2.0 


a = 


1 


> 12.2 






a = 


1/2 


> 11.7 






a = 


1/4 


> 2.8 



The best-fit cutoff energies are presented in Tabled Except for 
a - 1/4 (which actually implies very slow steepening of the 
electron spectrum rather than a cutoff), the synchrotron cutoff 
energies are high. In particular for the most realistic case with 
r = 2.0, the synchrotron cutoff energy exceeds 10 keV. 

3. Discussion 

A striking result of this study is the discovery of the inhomoge- 
neous, filamentary structure in the synchrotron X-ray emission 
of the NW rim. The bright filaments and hot spots embedded 
in the diffuse plateau emission, cannot be easily explained by 
the standard diffusive shock acceleration theory. More surpris- 
ingly, the X-ray spectrum does not show noticeable variations 
over the entire rim, in spite of the strong gradient of the surface 
brightness distribution. 

The derived hard power-law spectra of synchrotron X-rays 
imply that the electron acceleration continues effectively to 
> 100 TeV, unless we assume very strong magnetic field ex- 
ceeding 100 pG. Even so, independent of the strength of mag- 
netic field, the standard shock acceleration theory predicts a 
cutoff in the spectrum of synchrotron radiation around 1 keV 
or below, and correspondingly - steep X-ray spectra. 

In the framework of diffusive shock acceleration model, 
the synchrotron cutoff energy is set by the condition accelera- 
tion rate — synchrotron loss rate. The synchrotron cooling time 
of an electron of energy E e in the ambient magnetic field of 
strength B is 

fsvnch = 1.25 X 10 3 ( — — — ) [— — ) yr . (4) 
y \ 100 TeV/ \l0pG) y 

The diffusive shock acceleration time (see e.g. Malkov and 
Drurv 120011 can be written, with accuracy of about ±50 %, in 
a simple form f acc = 102D/V 2 , where T> is the diffusion co- 
efficient in the upstream region, and V is the upstream veloc- 
ity into the shock. The diffusion coefficient is generally highly 
unknown parameter, however, if one requires acceleration of 
particles to highest energies allowing synchrotron radiation up 
to X-ray energies, we must assume that the diffusion proceeds 
in the so-called Bohm diffusion regime, therefore it is conve- 
nient to parametrize the diffusion coefficient in terms of the 



Bohm diffusion coefficient Db = T]r g c/3, where r g = E e /eB 
and T] > 1 are the so-called gyroradius and gyrofactor, respec- 
tively. The value 77 = 1 implies the smallest possible value of 
the diffusion coefficient, and correspondingly the shortest pos- 
sible acceleration time, 

^ 2 - 7xl0 1i^)(i^) \moL^) 2 ^- {5) 

From Eqs. (0} and (|3J we find the maximum energy of acceler- 
ated electrons, 

^^Hl^rt^OOO^)^ 2 ^- (6) 

Now, substituting E m from Eq. (|6j into Eq. Q, we find that the 
synchrotron cutoff energy depends only on the shock speed and 
the gyrofactor, 

e ~ 2 ( T ) rf l keV . (7) 

boOOkms- 1 / ' V ; 

This simple treatment does not tell us much how far the 
electron spectrum can be extended beyond E m , and correspond- 
ingly how looks like the synchrotron spectrum in the cutoff 
region around eq. The assumption that the electron acceler- 
ation spectrum is described by the "canonical" Eq. is in- 
deed a rough approximation. In fact, the highest energy elec- 
trons from the cutoff region may have more complicated dis- 
tributions, which, for example, may contain a pronounced pile- 
up just proceeding the cutoff (e.g. Melrose & Crouch 1997 1. 
However, this may only shift the position of the synchrotron 
cutoff to lower energies (Aharonian 2000 1. Finally, we should 
note that the position of the synchrotron cutoff given by Eq. (0 
should be considered as an absolute upper limit because it does 
not take into account other possible losses. 

The results in the previous section indicate that the syn- 
chrotron cutoff in the radiation components observed from the 
NW rim, most likely are located above 10 keV. This requires 
V > 5000 J] l/1 km s" 1 which for any reasonable shock speed 
in RX J1713. 7-3946 is quite a tough condition and hardly 
could be fulfilled even when the acceleration proceeds in the 
extreme Bohm diffusion regime (77 = 1). A possible solu- 
tion to this difficulty could be found if one assumes that the 
electron acceleration and radiation regions are effectively sep- 
arated. Indeed, assuming that the acceleration takes place in 
sites of low-magnetic field, then the electrons quickly escape 
the acceleration region, enter into regions with significantly en- 
hanced magnetic field, and there produce the bulk of the ob- 
served X-ray flux with a synchrotron cutoff well above 1 keV. 
Otherwise, one should invoke faster (yet unknown) electron 
acceleration than the standard shock-acceleration model pro- 
vides, in order to account for the observed hard synchrotron X- 
ray spectra. Another possibility would be if the electrons are of 
secondary (7r ± -decay) origin produced at interactions of accel- 
erated protons and ions with the ambient gas. This hypothesis 
requires very strong magnetic field and acceleration of protons 
to energies > 10 15 eV. Such a model recently was suggested by 
Bell & Lucek (2001 1. This hypothesis implies very hard spec- 
trum of 7r°-decay TeV y-rays with energy flux comparable to 
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the energy flux of synchrotron X-ray emission. Therefore, the 
"hadronic" origin of electrons would be turned down, if the 
TeV y-ray spectrum indeed breaks at sub-TeV energies as has 
been claimed by Enomoto et al. ( 2002 1. 

Formally, the X-ray brightness distribution seen in Fig. @] 
could be result of different depths of the line of sight, i.e. as- 
suming that the bright filaments are sheets seen at "edge on", 
and the plateau emissions are those seen at "face on". Below we 
adopt a more likely, in our view, interpretation that the observed 
brightness distribution is a result of highly inhomogeneous pro- 
duction and distribution of relativistic electrons and magnetic 
fields in the NW rim. At first glance, homogeneous production 
of multi-TeV electrons over the NW rim cannot be ruled out, 
because in this case the filamentary structure might be referred 
to a local enhancement of the magnetic field. However, in the 
(most likely) regime when the synchrotron cooling time is less 
than the source age, the synchrotron X-ray luminosity does not 
depend on the strength of magnetic field. In this regime an equi- 
librium is quickly established between the electron injection 
and synchrotron losses. This implies that the synchrotron ra- 
diation saturates to the maximum possible rate determined by 
the electron injection rate, L sync h W e . Therefore, below we 
assume that the filamentary structure is associated, first of all, 
with inhomogeneous spatial injection of electrons. In general, 
since the propagation speed in different parts of the rim could 
be very different depending on the magnetic field and the devel- 
oped turbulence, we should expect significant spatial variation 
of the density of particles even in those cases when they are in- 
jected into the rim almost homogeneously. The acceleration of 
electrons in selective regions of the rim would make the X-ray 
synchrotron image of the rim even more irregular and clumpy. 

In this context, it is challenging to identify the X-ray fil- 
aments as sites of effective electron acceleration to TeV ener- 
gies. It is likely that the filaments are bright not only due to en- 
hanced magnetic fields, but also due to local concentrations of 
the multi-TeV electrons being directly accelerated there. To our 
knowledge, the filamentary structure is a feature which has not 
been discussed (dismissed ?) in the relatively well developed 
concept of diffusive shock acceleration in SNRs, although per- 
haps this could be attributed to the inhomogeneity of the ambi- 
ent medium. Also, it is interesting to note in this regard, that an 
appearance of "hot spots" or "discharge zones" has been quali- 
tatively predicted by Malkov ( 1997) who related such features 
to the drop of the critical injection caused by changes in the 
flow structure. Thus the very fact of existence of nonthermal 
filamentary structures of the shell of this SNR should not be 
interpreted as contradiction to the diffusive shock acceleration 
model. The unusually small spatial scales of these structures 
perhaps also could be accommodated by this theory assum- 
ing very large ambient magnetic field significantly exceeding 
10 fiG. However, as discussed above, any assumption of the 
magnetic field cannot explain the position of the synchrotron 
cutoffs beyond 10 keV, unless we assume more effective ac- 
celeration mode (e.g. Jokipii 1987 1. The question of the nature 
of the mechanism of particle acceleration is beyond the frame- 
work of this paper, therefore below we will adopt a priory an 
acceleration spectrum extending well beyond 100 TeV at the 
presence of relatively large field - a condition dictated directly 



by the X-ray observations - and explore the spectral and spatial 
features of synchrotron radiation associated with the propaga- 
tion effects of electrons in the rim. 

The small lateral dimensions of the compact structures in 
the filaments allow effective escape of accelerated electrons, 
thus the main part of electrons would eventually reside outside 
their production sites. Therefore, it is not a surprise that the 
total X-ray flux from the NW rim is contributed mainly by the 
plateau region, although the filaments look much brighter. 

Given the fast synchrotron cooling of X-radiating electrons 
compared to the source age and the compactness of the ac- 
celerators, any appropriate treatment of the problem should 
be time-dependent and should include propagation effects. 
More specifically, for the synchrotron cooling time r syn « 
760(e/5keV)" 1/2 (B/10/iG)' 3/2 yr, the spectral steepening of 
X-ray emitting electrons becomes unavoidable even for the 
minimum possible age of the source, To ~ 1000 yr. The 
energy-dependent escape on timescales r^f ^ AR 2 /2Ds - 
1500 dlT]-\E/ 10 TeV) _1 (B/10/iG) yr, is another important 
factor for modification of the electron spectrum in thin fila- 
ments. Here AR = 0.58 d(, pc is the typical thickness of fila- 
ments, and Db(E) = rjr g c/3 is the Bohm diffusion coefficient 
with the gyrofactor as a free parameter. In addition to the dif- 
fusive escape, we should take into account also the convective 
escape that occurs on timescales r con ~ pAR/V where p is the 
shock-compression ratio. The convection is an important effect 
for delivery of low-energy particles from filaments to the sur- 
rounding plateau region. 

To simplify the problem, and to clarify the basic relations, 
here we adopt a "two-zone" model applied earlier to SN 1006 
(Aharonian & Atoyan 1999 1. This model takes into account 
the effects related to the diffusive and convective escape of 
electrons from one homogeneous region (zone 1 : acceleration 
site, i.e. filaments) to another (zone 2: the rest of the rim, i.e. 
plateau, where electrons coming from zone 1 are accumulated 
without escape) with essentially different physical parameters 
such as magnetic field, diffusion coefficient, etc., and gives 
time-dependent solutions for energy distributions of electrons 
in both zones. In the Chandra image of the NW rim, the clas- 
sification of the filament/hot-spot (zone 1) and plateau (zone 2) 
regions is based on the surface brightness. The overall (i.e. in- 
tegrated over the volumes) fluxes in the 2-10 keV band are esti- 
mated to be l.lxlO" 11 ergcm~ 2 s and4.2xl0~ n ergcirT 2 s _I 
for zone 1 and 2, respectively. Figure|5]shows the X-ray spectra 
from zone 1 and 2 characterized by a power-law function with a 
photon index F = 2.2. Below, we try to reproduce the observed 
fluxes and spectra of synchrotron X-rays from the filament and 
plateau, simultaneously, in the framework of the "two-zone" 
model, in order to testify the "filament/hot-spot = accelerator" 
scenario. 

For the case of an old (tq = 10000 yr) and distant (d = 6 
kpc) SNR, we found that a satisfactory fit is possible assum- 
ing strong and similar magnetic fields in the filaments and 
plateau, = B p i a = 50 pG. Since the strong magnetic field 
implies fast synchrotron cooling, in order to get large X-ray 
fluxes from the plateau we should require an adequately fast 
escape of electrons from the filaments, namely r<jif < r syn , 
or 77 > 16 (Bfii/50 pG) 3 dl- In Fig.|9^ we assume 77 = 83 and 
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Tcon = 1000 yr. The X-ray slope of the filament is explained 
by synchrotron radiation from the diffusive-escape-steepened 
distribution of high-energy electrons, while that of the diffuse 
plateau - from the synchrotron-loss-steepened distribution. The 
derived parameters lead to quite similar X-ray spectra for two 
zones. The total energy of electrons accumulated in the plateau, 
Wpi a = 4.2 x 10 47 erg, appears one order of magnitude larger 
than that in the filaments. On the other hand, if the SNR is 
much younger (to = 1000 yr) and closer (d — 2 kpc), a good 
fit is possible assuming that the magnetic field is stronger in 
the filaments, namely = 20 pG, but weaker in the plateau, 
fipia = 6 pG. Since in this case the diffusive escape is very fast 
(Tdif °c AR 2 oc d 2 ), we must assume Bohm diffusion regime, 
r] = 1, in order to keep electrons within the filaments. We also 
assume r con = 500 yr. Due to the smaller magnetic field in the 
plateau region and the small To/r con ratio, the radio emissivity 
of the plateau is considerably suppressed. The radio flux com- 
ing from the NW rim region plotted in Fig.|9]should not exclude 
this model because the population of radio-emitting electrons 
and that of X-ray-emitting electrons could be very different; the 
high quality data of the radio fluxes would be needed prior to 
definitive conclusions. 

Modeling of the observed synchrotron X-ray emission is 
crucial for predictions of the spatial and spectral distributions 
of the associated IC y-rays produced by the same electrons up- 
scattering 2.7 K CMBR. In Fig.|9]we present the IC spectra 
from the filaments and the plateau. In both cases, we found that 
it is quite difficult to reproduce the flux and spectral shape of 
TeV y-rays reported by the CANGAROO collaboration, in the 
synchrotron-IC model. 

It should be noticed that the time-dependent treatment of 
the problem requires very effective acceleration of electrons. 
The exponential cutoffs in the acceleration spectrum are as- 
sumed at 125 TeV (case a) and 200 TeV (case b) in order to 
fit the X-ray data. As it follows from Eq. (jfji these assumptions 
hardly could be accommodated by the standard diffusive shock 
acceleration model. 

Finally, we would like to make a short comment concern- 
ing the possible association of the X-ray void with one of the 
three molecular clouds (cloud C, see Fig.^ reported by Slane 
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et al. ( fT999t in the proximity of RX J1713.7-3946. The den- 
sity of this cloud is estimated from CO observations (Bronfman 
et al. 119891 . n ~ 400 d b l cm 3 . If cloud C has physical asso- 
ciation with RX J1713. 7-3946, it could be an interesting site 
for the production of very high energy y-rays through inter- 
actions of electrons and protons with the ambient dense gas. 
In order to explain the bulk of the reported TeV luminosity 
(Enomoto et al. 2002), L y ~ 4 x 10 35 d\ erg s" 1 , by w°-decay 
y-rays, the total energy of protons in this cloud should be close 
to 10 50 iig erg. This exceeds, by 2 or 3 orders of magnitude, 
the total energy released in relativistic electrons, assuming that 
the particle accelerator has been operating during the age of 
the source: W e = r L e 8.8 x 10 47 and 5 x 10 47 erg for 
d = 6 and 2 kpc, respectively (see Fig. [9). Thus, the hypothesis 
of the 7r°-decay origin of the reported TeV flux would require 
a very large proton-to-electron acceleration ratio, significantly 
exceeding the canonical "100:1" value, given that only a frac- 
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Fig. 9. Multiwavelength synchrotron (solid lines) and IC 
(dashed lines) spectra of the filaments and the plateau re- 
gions calculated within the "two-zone" model, assuming that 
the electron acceleration takes place in the filaments. The 
ATCA flux of 4 Jy at 1.36 GHz (squares) is from Ellison 
et al. ( I2U01I . The X-ray spectra from the filaments and the 
plateau region are from this work, and the TeV y-ray flux is 
from the CANGAROO observations (Enomoto et al. 2002 1. 
The EGRET upper limit corresponds to the flux of 3EG 
J1714-3857 (Butt et al. 1200 1> . The following parameter sets 
have been used in calculations: a) to = 10000 yr and d = 6 
kpc. r con = 1000 yr. Acceleration spectrum of electrons with 
the spectral index s — 1 .84, exponential cutoff at Eq — 125 TeV, 
and acceleration rate L = 2.8 x 10 36 erg s _I . fi^i = B p i a = 50 pG 
and Tj = 83. b) t = 1000 yr and d = 2 kpc, r con = 500 
yr, s = 1.95, E = 200 TeV, and L = 1.6 x 10 37 erg s -1 . 
Bjii = 20 pG, Bpi a = 6 pG and r\ = 1 . 



tion of the protons accelerated by the supernova shock could be 
captured by cloud C. 

4. Conclusions 

The Chandra image has revealed that the synchrotron X-ray 
emission from the northwestern rim of SNR RX J1713. 7-3946 
has remarkable fine-structure: the complex network of syn- 
chrotron X-ray filaments surrounded by fainter diffuse plateau, 
and a dark region with a circular shape. By examining individ- 
ual spectra, we found that despite significant brightness vari- 
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ations, the spectral shapes of the X-ray spectra everywhere in 
this region are more or less similar, being well fitted with a 
power-law model of photon index r - 2.3. 

The observed hard power law requires rather high syn- 
chrotron cutoff frequency (energy), set by the condition "ac- 
celeration rate = synchrotron loss rate", which is most likely 
to be > 10 keV taking account of the effects of spectral steep- 
ening due to synchrotron losses. We need unreasonably high 
shock speed exceeding 5000 km s _1 to explain such a high cut- 
off energy within the standard formalism of the diffusive shock 
acceleration model. A possible solution to this difficulty could 
be obtained if one assumes low magnetic field in the accel- 
eration sites, and high magnetic field in the emission regions. 
Otherwise we should invoke faster, as yet unknown, electron 
acceleration mechanism. 

We testify the scenario that the filaments/hot-spots are 
the acceleration sites of bulk of multi-TeV electrons without 
specifying acceleration mechanism, by adopting a "two-zone" 
model. Then we found that our time-dependent treatment in- 
cluding the effects of fast energy losses and diffusive escape 
of electrons is capable of accounting for the fluxes and spec- 
tral shapes of X-ray emissions from both the filament and 
plateau regions. We note that the TeV spectrum reported by the 
CANGAROO collaboration cannot be readily explained by the 
IC scattering of the CMBR photons off X-ray-emitting multi- 
TeV electrons. Higher quality data from the radio, X-ray and 
y-ray bands should help us to draw more definitive conclusions 
concerning the nature of acceleration mechanisms and the ori- 
gin of the observed TeV y-rays. 

Acknowledgements. We thank the referee, B. Aschenbach for very 
useful comments which helped to improve the paper, as well as 
G. Rowell and Y. Butt for helpful comments. YU. is supported by 
the Research Fellowships of the Japan Society for the Promotion of 
Science for Young Scientists. 



Reimer, O., & Pohl, M. 2002, A&A, 390, L43 
Slane, P., Gaensler, B. M„ Dame, T. M„ et al. 1999, ApJ, 525, 357 
Uchiyama, Y, Takahashi, T, & Aharonian, F. A. 2002, PASJ, 54, L73 
Weisskopf, M. C, O'Dell, S. L., & van Speybroek, L. P. 1996, 
Proc. SPIE, 2805, 2 



References 

Aharonian, F. A., & Atoyan, A. M. 1999, A&A, 351, 330 
Aharonian, F. A. 2000, New Astronomy, 5, 377 
Aharonian, F. A. 2001, Space Sci. Rev., 99, 187 
Bell, A. R., & Lucek, S. G. 2001, MNRAS, 321, 433 
Bronfman, L., Alvarez, H., Cohen, R. S., & Thaddeus, P. 1989, ApJS, 
71,481 

Butt, Y. M., Torres, D. F, Combi, J. A., Dame, T, & Romero, G. E. 

2001, ApJ, 562, L167 

Butt, Y. M., Torres, D. F, Romero, G. E., Dame, T, & Combi, J. A. 

2002, Nature, 418, 499 

Ellison, D. C, Slane, P., & Gaensler, B. M. 2001, ApJ, 563, 191 
Enomoto, R., Tanimori, T, Naito, T, et al. 2002, Nature, 416, 823 
Freeman, P. E., Kashyap, V., Ronsner, R., & Lamb, D. Q. 2002, ApJS, 

138, 185 
Jokipii, J. R. 1987, ApJ, 313, 842 

Koyama, K., Kinugasa, K., Matsuzaki, K., et al. 1997, PASJ, 49, L7 
Malkov, M. A. 1997, ApJ, 491, 584 

Malkov, M. A., & Drury, L. O'C. 2001, Rep. Prog. Phys., 64, 429 
Melrose, D., & Crouch, A. 1997, Proc. Astron. Soc. Aust., 14, 251 
Pfeffermann, E., & Aschenbach, B. 1996, in Rontgenstrahlung from 

the Universe, ed H. U. Zimmermann, J. Triimper, & H. Yorke 

(MPE Rep. 263) 267 



This figure "figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0209217v2 



This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0209217v2 



This figure "fig4a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0209217v2 



This figure "fig4b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0209217v2 



This figure "fig5.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0209217v2 



This figure "fig6.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0209217v2 



